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Abstract 

We analyze the behavior of a quantum system described by a one-dimensional asymmetric poten- 
tial consisting of a step plus a harmonic barrier. We solve the eigenvalue equation by the integral 
representation method, which allows us to classify the independent solutions as equivalence classes 
of homotopic paths in the complex plane. We then consider the propagation of a wave packet 
reflected by the harmonic barrier and obtain an expression for the interaction time as a function 
of the peak energy. For high energies we recover the classical half-period limit. 
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I. INTRODUCTION 



The harmonic oscillator plays a central role in physics because it is exactly solvable and 
provides a simple model for a host of physical phenomena. Besides the simplest case of the 
one-dimensional free oscillator, modifications of the harmonic oscillator are of interest. A 
typical variant is to confine the oscillator in a box. This system has been investigated in 
one dimension^- and for arbitrary dimensions.- The purpose of these investigations is to 
calculate the corrections to the energy levels caused by the presence of infinite barriers at a 
finite distance. This analysis has also been done for spherically symmetric potentials, such 
as the hydrogen atom.- In addition, the behavior of a wave packet propagating in a generic 
power-law one-dimensional potential well has been considered in terms of its "collapse and 
revival," namely of its scattering over the well and its subsequent reforming.- 1 ^ The truncated 
harmonic oscillator has been used to study how the presence of discrete levels in the energy 
spectrum affects tunneling through such a well. 9 

In the present paper we study the bound states and the propagation of a wave packet in a 
one-dimensional potential consisting of a half-space harmonic oscillator plus a step. In Sec. [XT] 
we solve the Hamiltonian eigenvalue equation using the integral representation method, 
which allows us to classify the independent solutions as equivalence classes of homotopic 
paths in the complex plane. In Sec. IIHI we calculate the energy of the bound states and 
compare them with the standard harmonic oscillator and the half-space harmonic oscillator 
with an infinite barrier. In Sec. IIVI we study the properties of the propagation of a wave 
packet which, coming from infinite distance, is reflected by the harmonic barrier. Our 
analysis is based mainly on the investigation of the delay time in the reflection, which can 
be interpreted as the duration of the interaction with the harmonic barrier. We express the 
interaction time as a function of the peak energy and study its asymptotic behavior. We 
show that in the high energy limit the delay approaches the classical value, namely the half 
period of the harmonic oscillator. Finally, we comment on how this behavior changes in the 
presence of a stronger or weaker confinement. 

The problem that we address has been studied by Mei and Lee^ to test the adequacy of 
a perturbation scheme on an exactly solvable model. Our analysis, which has a different 
purpose, has the advantage of being based on the integral representation method which can 
be applied to a wider class of problems. In addition, we do not confine ourselves to the study 
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of the bound states, but also investigate the motion of continuous spectrum wave packets. 



II. THE STEP-HARMONIC POTENTIAL 



Consider a particle subject to the one- dimensional potential: 



U (x > 0) 

U{x) =\ x (1) 

-KX 2 (x < 0), 



where k and Uq are real positive constants (see Fig. [p. 



FIG. 1: The step-harmonic potential. 



The proper and improper eigenfunctions of the Hamiltonian operator are ordinary solu- 
tions of the eigenvalue equation outside of the discontinuity of the potential. Such solutions 
must be continuous together with their first derivatives across the sing ularilyjiBill We will 
solve the eigenvalue equation for x > and x < 0. 

For x > the eigenvalue equation is 

h 2 d 2 u(x) TT , . _ , s 

where m is the particle mass and E is the energy eigenvalue. The general solution of Eq. ([2]) 
has the form: 

{ Ae ikx + Be -ikx ^ E > Uq } 
(3) 
A'e kx + B'e~ kx (0 < E < U ) 

where Hk = yj2m\E — U$\. 

We must choose A' = 0. Otherwise, for < E < Uq, u(x) would diverge exponen- 
tially for x — > +oo and therefore it would neither belong to L 2 (IR) (i.e. the space of 
the square summable functions over R), nor satisfy the eigenpacket condition for improper 
eigenf unct ions .— ^ 



For x < the eigenvalue equation is 

d 2 u(y) 

dy 2 

where 



+ (e-y 2 )u(y) = 0, (4) 



rriK 2E m Ik 

We set u(y) = F(y) exp (—y 2 /2) and obtain from Eq. (jlj) the following equation for F(y) 

F"(y)-2yF'(y) + (e-l)F(y) = 0, (6) 

which is the Hermite equation.— ^ 

The solutions of Eq. ([6]) are entire functions and can be found by the method of integration 
by series. We prefer to employ the integral representation method.— We start by looking 
for solutions of Eq. that have the form 

F(y)= [ dtf(t)e~ t2+2ty , (7) 



where 7 is a path in the complex plane C and / is a suitable function which is holomorphic 
in a region which contains the graph of 7. We substitute Eq. (J7|) into Eq. (Q and obtain 

J dt [At 2 + (e - \)}f(t)e- t2+2ty - 2 J dt (J^e 2ty ) tf(t)e~ t2 = 0. (8) 

Equation ([81) , after integration by parts of the second integral, can be written as 

-2t/(t)e-* 2+ H + [dt [(e + l)/(t) + 2tf'(t)]e' t2+2ty = 0. (9) 

From Eq. (jSJ) it follows that Eq. ((Tj) is a solution of Eq. Q) if 

tf{t)e- t2+2ty ] =0 and /(t) = t'^ . (10) 

Therefore, we can write a solution of Eq. ([6]) in the form 

F^\y) = [ dtr^e- t2+2ty , (11) 

where 7 must be chosen according to the first condition in Eq. (|T0|) and such that the integral 
in Eq. (TTTT) is well defined. The classification of the appropriate 7's allows us to classify all 
solutions of Eq. (JH]). 

The integrand in Eq. ( TTTT) is singular at t = 0. For e = 2n + 1 (n = 0, 1, . . . ), the point 
t = is a pole of order n + 1, otherwise it is a branch point. 

In the following we distinguish two classes of paths, which correspond to two linearly 
independent solutions of Eq. <M). 



FIG. 2: Possible paths for e = 2n + 1 (n = 0, 1, . . . ). 



A. The case e = 2n + 1 



In this case we can rewrite Eq. (j7J) as 



^ 7) (Z/) = / 




+2*2/ 



(12) 



where the integrand is holomorphic on C but the origin. Possible choices of 7 for which 
the contour condition in Eq. (flOl) holds are shown in Fig. [21 Ti and T3 have a real part 
that goes to infinity, T 2 is a closed path circling the origin, and T 4 is a closed path that 
does not contain the origin. By virtue of Cauchy's theorem, F^f 1 = 0, and because the 
paths can be deformed so that Ti + T3 = T2, the other three solutions satisfy the relation 
En + En = Fn \ where Fn^ is the solution corresponding to the path Tj (j = 1,2,3,4). 
Then we have, as expected, two linearly independent solutions for Eq. (JB]). 
As an exercise, we show that the solution corresponding to T 2 , namely 



corresponds to the Hermite polynomial of order n. By completing the square in the integrand 
of Eq. (EE3), we find 




(13) 




(14) 



We take advantage of Cauchy's formula and rewrite Eq. (I14j) as 



F< 2 'M = ^(-l)"^(e-»') = ^ 
n\ av 11 n! 



H n (y), 



(15) 



where H n (y) is the Hermite polynomial of order n.— 



Im(t) 




FIG. 3: Possible paths for e ^ 2n + 1. Note that in this case T2 cannot be closed at infinity. 
B. The case e / 2n + 1 



In the generic case e G K, e ^ 2n+ 1, we rewrite Eq. ( II ID as 



-t 2 +2ty 

F^(y)= I M—^, (16) 
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where = (e + l)/2. If f3 is not a positive integer, t = is a branch point for the multivalued 

function t@ . In this case we must cut the complex plane, for example along the positive real 

axis. In the latter case, the classes of possible paths are depicted in Fig. |3j 

In the following we show that the solutions corresponding to Ti and T 3 (iv and F^) 

diverge as e y2 for y = ±00 and therefore the corresponding eigenfunction u(y) cannot be 

(2) 

either proper or improper. We are thus left with the solution corresponding to T2 [iv ] which 
again diverges as e y2 for y — > +00. It also diverges for y — > —00, but the corresponding u(y) 
and its derivatives vanish more rapidly than any polynomial thanks to the presence of the 
exp(— y 2 /2) factor. 

According to standard results in the theory of integrals depending on a parameter, it is 
easy to show that F e (y) is an entire function. Furthermore, the derivatives of are 
obtained by differentiating with respect to y under the integral sign of Eq. ( TL6l) . Thus we 
obtain the relation 

We next address the asymptotic behavior of two independent solutions, for example, F e 
and F e (2) . 
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The Ti solution. We rewrite Eq. f|T6|) for the path I\ by introducing the variable z = t—y: 



_,2 

e 



F > 1>{ ^ e " k iz WTW (18) 

The branch point is now z = —y and the cut is shifted as well. We extract \y\ from the 
integral in Eq. ( fl8|V ^: 



ly P\ 



FP(V) = UmG 1 (y) 1 (19) 

where 

/ ri [sgn(y) + A] 



_,2 

e 



Gi(v)=/ ^ r _^ A , (20) 



An elementary calculation shows that 



e" 2 ' x/tt 



lim GAy) = lim / dz 7 — — r = (21) 

^±00 uy; ^±-7 ri [sgn(y) + ^<VY 

where we have taken the limit under the integral sign by virtue of the dominated convergence 
theorem^ (see Appendix [B]). The asymptotic behavior of F^\y) for y — > ±00 is therefore 

i^O/H-v 7 ^. (22) 



The corresponding eigenfunction u^(y) = Fe (y) exp (— y 2 /2) cannot be either proper or 
improper. 

(2) 

The T 2 solution. The solution corresponding to T 2 [Fe (y)\ has the form 

r p -t 2 +2ty 

Fp\y)= tot-p-, (23) 

where, as shown in Fig. |3l T 2 circles around the branch point in an anti-clockwise sense. 
This solution has different behavior for y — > +00 and y — > —00. 

Asymptotic behavior for y — > +00. By virtue of Cauchy's theorem we can deform r 2 to 
split the integral in Eq. ( )23l) into a sum of two integrals over the paths Ti and T^: 

e -t 2 +2ty 



f?\v)= dt ^?—- (24) 

Jriur 3 1 

We again introduce z = t — y and extract y 13 from the integral. We obtain 



F?\y) = ^G 2 {y\ (25) 
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Im(t) 



Si 



Re(t) 



FIG. 4: A standard trick in contour integration. 



where 



y> "'3 v y 

Thanks to the dominated convergence theorem we find 

lim G 2 (y) = -2ie-^V^sin(7r/3). 

y— >+oo 
(2) 

The asymptotic behavior of _F e (y) for y — > +oo is thus 



F^(y) ~ -2ie~^ /3 v ^Fsin(7r/3)^. 



(26) 



(27) 



(28) 



Note that, if e = 2n + 1, Eq. ( 1281) is incorrect because G2(y) — >• for y — )■ +00. We already 
know that, in this case, F ( - 2 \y) is as a polynomial of degree n. 

Asymptotic behavior for y — > —00. By taking advantage of Cauchy's theorem, we deform 
and split T 2 in the 3 sub-paths shown in Fig.Hl We choose for simplicity r = 1; from Eq. (123"]) 
we obtain 

poo ~t 2 +2ty 

F®{y) = I fi {y) - 2«f** sin(7r/3) jf dt — ^ , (29) 

where 

p2n 

Ip(y)=i £Q e i(l-l3)0 e -co S (20)-ism(2e)+2y cos 0+2iy sin0 _ ^q-j 

For y < and t > the following inequality holds 



t^e- t2+2ty < t~V 2 , 



(31) 



and, moreover, t _/3 exp(— t 2 ) is integrable in [l,+oo). Therefore, the integral on the right- 
hand side of Eq. (1291) vanishes for y — > —00 by virtue of the dominated convergence theorem. 
In contrast, for Ip(y) we have 



•2tt 

< / d6 





1(1-/3)0 -cos(20)-isin(20)+2j/ cos6+2iy sing 



i-2-k 

/ d^e- cos(2e) e 2?/cose , (32) 



so that < 27re 2 l y l 1 . Therefore, for y — > — oo, the absolute value of F { e (y) is domi- 

nated by 2ne 2 \ y \~ 1 , and 

u ll {y) = F^\y)e-y 2 ' 2 (33) 

is rapidly decreasing. Hence it is square summable on the positive real axis. Thus, it 
follows that for e 7^ In + 1 the full-space harmonic oscillator does not admit proper or 
improper eigenfunctions. More general theorems^- 2 allow one to obtain our results indirectly, 
for example, by studying the asymptotic behavior of the power series expansion of the 
solutions of Eq. Here we have adopted a more direct approach. 



III. EIGENFUNCTIONS AND ENERGY LEVELS 

From the results of Sec. [TT] we can write the energy eigenfunctions as 

{«V 
AFJax)e 2 ( x < 0) 
(34) 
Be~ kx (x > 0) 

{CFAax)e~^f~ (x < 0) 
(35) 
De ikx + Ee -ikx ( x>0 ), 

where we have dropped the superscript from F e . The integration constants A,...,E must 
be chosen such that u(x) and its first derivative are continuous at x = (the junction 
conditions) . 



A. The case E < U 

If the energy is smaller than the step height Uq, the junction conditions imply that 

B-F e (0)A = (36a) 
kB + aF^0)A = 0. (36b) 

The condition for the existence of a nontrivial solution is the vanishing of the system de- 
terminant. We define J{(3) = F e (0) [see Eq. (IA1B ]. In Appendix lAl we obtain the following 
expression for J(f3) [see Eq. fl A3[) ] : 

^)^r(i^). (37 , 
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The recurrence relation for the derivatives of F e , Eq. ( ITT)) , can be rewritten in terms of J ({3) 
as 

F e '(0) = 2J(/3-l). (38) 
Equation (138 j) implies that the junction conditions can be rewritten as 

-2aJ(/3- 1) = fcJ(/3), (39) 

or, equivalently, as 



i-2 



r 



1-/3 



2/ //3o-/5 



(40) 



2 

where /3q = Uo/(hw) + 1/2. By using the relation- 



we can rewrite Eq. (J40J) as 

r 



rwr(i -z) = ^— (4i) 

sm(7rz) 



/3 + 1 



cot(^)=-J^, (42) 



For a given value of /3q Eq. (142)) is an implicit relation determining the energy levels. The 
advantage of Eq. (142)) is that the singular behavior is contained in the cotangent function. 

For < Uq < hw/2 (1/2 < (3q < 1) the step is too small to allow for the existence of 
discrete energy levels. The first level (the ground state) appears for (3 = 1 at the value 
Eq = hw/2, which is the ground state of the full-space harmonic oscillator. For 1 < (3q < 3 
there is only one level, whose energy grows with increasing (3q starting from its minimum 
value hw/2. When (3 crosses the value 3 a second level appears at the energy Ei = bhw/2. 
By further increasing f3 (and hence Uq) the subsequent levels appear as Uo crosses the values 
Ek = hw(2k + l/2) (k e N), corresponding to the (fc + l)th even level of the oscillator. Thus, 
for a fixed value of /3q such that 2k + 1 < O < 2k + 3 there are exactly fc+1 bound states with 
energies E n (n — 0, 1, . . . , k) satisfying the inequalities hw(2n + 1/2) < E n < hw{2n + 3/2). 
An example with k = 1 is shown in Fig. [5j Each E n is a monotonically increasing function 
of Uq which asymptotically approaches the value hw(2n + 3/2) (the (n + l)th odd state of 
the oscillator) as U —¥ oo. 
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FIG. 5: The solid and the dashed lines represent respectively the left- and the right-hand side of 
Eq. (|42p . The intersections determine the energy levels. Here /3q = 4.5. 

The (unnormalized) eigenf unctions corresponding to the eigenvalue E n are 

F en {ax)e~^T ( x <0) 
u n (x) = { (43) 

J(f3 n )e~ k - X (x > 0), 



where the f3 n are the solutions of Eq. ( 1421) ; e n = 2(3 n — 1, Hk n = ^/2m([/ — E n ), and 
E n = hwe n /2. 

B. The case E > U 

The junction conditions on the eigenfunctions of Eq. (135]) are 

D + E — CF e (0) = (44a) 
ik(D — E) — CaF e '(0) = 0, (44b) 

implying the normalized (with respect to k) improper eigenfunctions are given by 

( «V 

1 U(l3)FJax)e~ 2 ( x < 0) 
Ue (x) = -^={ ^ (45) 
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where, as usual, 2/3 = e + 1, Hk = \/2m(E — Uq), 2E = hue, and 



n(/3)=2 j(/3) + i A /^-J(/3-l) 



C(/3) 



r 



1-/3 



/3-/V 



2-/3 



r 



1-/3 



/3-/V 



2-/3 



(46) 
(47) 



where we have used Eq. (I A 1 [) . Note that |C(/3)| = 1. As expected, the continuous part of 
the spectrum (E > U$) is simple. 



IV. REFLECTION AND DELAY 



To study reflection phenomenon, we consider the following superposition of continuous 
states: 

il){x,t) = j dkc{k)u <k) {x)e'h mt . (48) 
Jo 

From Eq. (T4"5"]) we have 



i/)(x,t) 



1 | J^dkc(k)U^(k))F e (ax)e—^-n E{k)t 



{x < 0) 



271 1 r o °° dA; c(k) [C(f3(k))e ikx + e~ ikx ] e" i_ ^* = ^ ref + (a > 0) 
We write ip in and ip Te f in the form: 

1pref(x,t) -■ 

where we have defined 



(49) 



1 

J /"+OC 

/2T 







dfc|c(A;)|e- i[M(fc) '- 7(fc)] , 
dk \ c (k)\e i[kx ^ k)t+s( - k)+ ' 1< - k) \ 



e im = mk)) and fi(fc) ^^P = ^ + M2 



(50) 
(51) 

(52) 



h h 2m 

If c(k) is sufficiently regular and non-vanishing only in a small neighborhood of k, then ip 
and ^ref represent wave packets that move according to the equations of motion^^ 

d7 



dn 

dk 



t + 



k=k 



dk 



hk , s p , . 
(t- t ) = --(*- to) 



fc=fc 



m 



m 



for the "incoming" wave packet, and 

•^ref 



dfi 


t - 

k=k 


d7 




d8 


_ _p 


dA 7 


dfc 


k=k 


dk 


k=k m 



m dS 

(t — to) ~ ~jT 

p dk 



k=k. 



(53) 



(54) 
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FIG. 6: Plots of the delay time r (in units of T/2, where T = w n/m is the period associated to 
the harmonic oscillator) versus the "energy" j3 of the incoming wave packet for six values of the 
step height A), (a) /3 = 1.5, (b) ft = 2, (c) /3 = 2.5, (d) /3 = 3.5, (e) /3 = 4, and (f) ft = 4.5. 

for the reflected "outgoing" one. 

The solution represents a particle of well defined momentum p = hk which approaches 
the origin from the right, interacts with the harmonic potential (at t = to), and is totally 
reflected. The phase shift results in a delay in the time the wave packet bounces back, which 
is caused by the interaction with the confining harmonic barrier. Because the phase shift 5 
depends only on k through (3, we can write the delay as 



where (3 = (3(k). 
We prove that 



1 d5_ 

u d/3 



(55) 



lim 5'ifi) 



lim S'((3) = 7r, 

+oo (3 e U(2A;,2A; + 1) 

fceN 

-oo fa E U (2k + 1,2k + 2) 

k&N 

/3 GN = 1,2,..., 



(56) 



(57) 
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where the prime denotes the derivative with respect to (3. To this purpose, note that by 
employing Eq. ( 14 ip 5'{0) can be cast in the form: 



HP) 



^v^^sin^vr) 



1 



/3-/3 



P 



+ 



2tt 



sin (/?7r) 



(P-Po)T(P/2) 
T (0/2 + 1/2) a/2 



sin 



/3vr 

T 1 



r(/3/2 + l/2) y^ 

r 03/2) 



cos^ 



(58) 



where $ is the Digamma function (that is, the logarithmic derivative of the Gamma 
function).— In Fig. E] we plot r versus (3 for different values of /3q. Note the resonances 
located at /3 ~ 3, 5, 7, 9, ... , corresponding to the formation of metastable states at the re- 
spective energies E ~ 5ftw/2, 9hu/2, 13fuj/2, 17hw/2, .... These states have lifetimes which 
decrease as the corresponding energies increase and move farther away from the threshold 
energy Uq. Conversely, as [7 increases, the lifetime of the resonance closest to the height of 
the step becomes progressively longer and then infinite when the resonance turns into the 
next bound state. This behavior is evident in Fig. [HI in which the first three plots correspond 
to values of fto for which there is only one bound state. In the successive three plots the 
resonance at (3 = 3 has disappeared, having turned into the second bound state. 

It is simple (using steepest descent or Stirling's formula, for example) to show that 



T{z + 1/2) 
T(z) 



1 + ( - 

z 



(59) 



for z ^> 1. By using one of the integral formulas for the Digamma function, 13 we can also 
show that 



lim 

z— >oo 



* (z) - * [Z + 



0. 



(60) 



Thanks to Eqs. ( |59j) and ( |60l) it is straightforward to derive Eqs. (1561) and ( 1571) . In particular 
Eq. (156|) implies that 

7T T 

(61) 



7T T 

lim t(8) = - = -. 



The wave packet undergoes half an oscillation during the interaction with the harmonic 
potential before being reflected, which results in a delay of half a period compared with 
the reflection from a perfect mirror (that is, when the confining barrier is an infinite wall). 
Thus, as expected, the high energy limit reproduces the classical behavior. 
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V. CONCLUSIONS 



The main features of the discrete part of the spectrum can be summarized as follows. 
For sufficiently small Uq (the height of the step) there is no discrete spectrum. When Uq 
increases and approaches the value hu/2 from below, there appears a resonance at energy 
Eq ~ hu/2. This resonance converts into a bound state when Uq reaches the value huj/2 
and the corresponding eigenfunction is proportional, at x < 0, to the eigenfunction of the 
ground state of the free harmonic oscillator and is flat otherwise. By further increasing 
the height of the step the ground state energy increases monotonically with Uq and, as 
Uq — > oo, approaches asymptotically from below the first odd level 3hco/2 of the full-space 
harmonic oscillator. A new discrete energy level appears at each energy Ef. = hu(2k + 1/2) 
(k e N) whenever Uq crosses the value E^. In the limit of infinite Uq (leading to the half- 
space oscillator), the energy levels become the odd levels of the oscillator itself, as expected 
from the symmetry of the problem. Loosely speaking, the levels are born as "even" and, 
upon increasing the height of the step, end up as "odd" (see Fig. [5]). This behavior is not 
peculiar to the step problem associated with the harmonic oscillator, but is typical of the 
corresponding step variant of every symmetric confining potential. 

The continuous spectrum is simple and extends from Uq to oo. A wave packet coming 
from infinity collides with the confining harmonic branch and is thereby entirely reflected. 
The interaction with the potential results in a delay of the reflected packet which, as is 
well known for problems of this kind, is proportional to the derivative of the phase shift 
of the plane wave component evaluated at = f3(k). This delay can be interpreted as the 
interaction time with the harmonic barrier. When the confining part of the potential is 
infinite {U(x) = +oo at x < 0) the delay vanishes, and the reflection on a perfect mirror is 
instantaneous. The desirable feature of our example is that we can derive an exact analytic 
expression for the delay [Eq. fl58l) ] as a function of the step height and of the peak energy 
of the incoming packet. 

Although these characteristics are typical of the step variants of all symmetric confining 
potentials, the harmonic oscillator potential is "more equal" than the others. It is the only 
analytic (except possibly at x = 0), convex or concave locally bounded symmetric and 
confining potential that gives rise to classical isochronous oscillations and thus to evenly 
spaced energy levels.— 1 ^ In our step variant of the problem we recover both these features 
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in the limit Uq — > oo, when the potential reduces to the half-space harmonic oscillator. Thus, 
it is possible that the harmonic potential is the only confining barrier that displays a constant 
nonvanishing interaction time in the limit of high energies. For steeper barriers we expect the 
interaction time r to vanish at high energies, while for milder potentials we expect the delay 
to become infinite in this limit, in accordance with the corresponding classical behavior. 
Similarly, we expect that, as Uq — > oo, the spacing between two neighboring discrete levels 
approaches infinity in the former case and zero in the latter. 

In a forthcoming paper we corroborate this conjecture by analyzing two examples using 
the integral representation method. This method can be employed to analyze a wide class 
of "step-something" potentials in which the harmonic part is replaced by another type of 
barrier. We encourage readers to investigate, for example, the step-linear (si) and the step- 
exponential (se) potentials 



where M, a, and Uq are real positive constants. 
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Appendix A: Calculation of </(/?) 

In this appendix we prove Eq. ( 13"T|) . According to the definition given in Sec. [HI] we have 



where 2(3 = e + 1. Assume f3 G C. For \f3\ < R, the integrand function in Eq. (IA1|) is 
bi-continuous and holomorphic with respect to in any compact disc. Furthermore, its 
absolute value is bounded by a summable positive function: |e _ * 2 t _/3 | < e~ 2 ^ t2 ^\t\ 2R . These 




(62) 



(63) 




(Al) 
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FIG. 7: The path T2 and its transformed one V 2 after the change of variable u = t 2 . When we take 
into account the presence of the cut, we must choose < arg(u) < 2ir on the first sheet (dashed 
line), and 2ir < arg(n) < 4n on the second one (solid line), that is, we choose the positive square 
root on the first sheet, and the negative one on the second. 

properties imply that the integral in Eq. flAlj) is uniformly convergent. Therefore, J(/3) is 
an entire function. If we change the variable in Eq. (lAlj) to u = t 2 , we must cut the complex 
plane and define t = sfu on its complete Riemann surface, which is composed of two sheets. 
The new path T' 2 is shown in Fig. [7J 

On both sheets of the w-plane, there is an integral along a straight line and an integral on 
a semi-loop. It is easy to show that the integral on the semi-loop vanishes when shrunk to 
a point, provided that G (— 00, 1). The integrals along the straight lines can be deformed, 
using Cauchy's theorem, to become integrals along the positive real axis. We thus have 



JO?) 



which we can write as 



If 00 1 i4n 

due u u 2 + - e 



2 



3 _ _£±1 
due u u 2 , 



ic. 



(A2) 



(A3) 



Because the poles of the Gamma function in Eq. ( 1A3I) are cancelled by the zeroes of the 
sine, the right-hand side is an entire function so that Eq. 0A30 holds on the entire complex 
plane by analytic continuation. 



Appendix B: Taking limits under the integral sign 

For convenience we give here a very useful elementary theorem which we have used in 
the paper: 

Theorem 1: The dominated convergence theorem. Let fk '■ K — > K be summable 
on an interval /, that is, fj fk < 00, VA; 6 N. Moreover, let fk converge almost everywhere 
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to a function : R — > R. Suppose that there exists a positive J-summable function 
g : R — > R + that dominates every f k (that is, < g(x)Vk e N). It follows that /oo is 

J-summable and that we can take the limit under the integral, that is, 

lim / dxf k (x) = dx lim f k (x) = / dxf^x). (Bl) 

k^oojj Jj k-^oo Jj 

It is noteworthy that, in our case, we do not need to invoke Lebesgue integration and the 
dominated convergence theorem, and our calculations are based on Riemannian integration. 
Even though Riemann integration theory lacks theorems regulating the interchange between 
limit and integration operations, the following theorem is sufficient for our purposes:— 

Theorem 2. Let f n be a sequence of functions defined on [a, oo) and Riemann integrable 
on [a, b] for all b > a. Assume that (i) f n {x) —> foo{x) almost everywhere in [o, oo), foo being 
Riemann integrable on every finite interval, (ii) There exists a positive function g defined 
on [a, oo) such that f°°g is convergent and |/ n (^)| < g(x) for all n. Then f n — > f a °° f^. 
Note that, in this theorem, the integrability of the limiting function is part of the hypothesis, 
whereas in the dominated convergence theorem it is a consequence of the theorem itself. 
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